Comparison of 32-site exact diagonalization results and 
ARPES spectral functions for the AFM insulator Sr 2 Cu0 2 Cl2 
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We explore the success of various versions of the one-band 
t — J model in explaining the full spectral functions found 
in angle-resolved photoemission spectra for the prototypical, 
quasi two-dimensional, tetragonal, antiferromagnetic insula- 
tor Sr2Cu02Ci2. After presenting arguments justifying our 
extraction of A(k,u>) from the experimental data, we rely 
on exact-diagonalization results from studies of a square 32- 
site lattice, the largest cluster for which such information is 
presently available, to perform this comparison. Our work 
leads us to believe that (i) a one-band model that includes 
hopping out to third-nearest neighbours, as well three-site, 
spin-dependent hopping, can indeed explain not only the dis- 
persion relation, but also the quasiparticle lifetimes — only in 
the neighbourhood of k = (■§, 0) do we find disagreement; (ii) 
an energy-dependent broadening function, T(E) =Fq+ A E, 
is important in accounting for the incoherent contributions to 
the spectral functions. 



I. INTRODUCTION 

The study of strongly correlated electronic systems in 
two dimensions (2D) has been very active, particularly 
since the discovery of the high T c superconductors. How- 
ever, while there have been great improvements in the 
reliability and scope of experimental measurements, this 
has not led to agreement as to the correct description of 
either high temperature superconductivity, in particular, 
or two-dimensional correlated systems, in general. Quite 
the opposite is in fact apparent — certain experimental 
results are consistent with one phenomenology, but not 
all experiments seem to support the same physical con- 
cepts, and as of yet no one theory is capable of explaining 
all the data. One avenue that may alleviate this confu- 
sion and allow for a consensus to develop follows from 
the study of ideal experimental systems. By this we im- 
ply certain compounds are amenable, for various physical 
as well as circumstantial reasons, to a complete charac- 
terization. Further, these compounds are such that a 
relatively simple theoretical model, believed to be appro- 



priate for that material, can be defined, and with suitable 
ingenuity we can expect to carry out a fully quantitative 
comparison of this model and the experimental data. 

A potential paradigm of 2D antiferromagnetic (AFM) 
insulators, a category which includes the parent com- 
pounds of high T c superconductors, is the Sr2Cu02Cl2 
compound [I]]. The features that make this compound 
ideal are: (i) it remains tetragonal down to 10 K; (ii) it 
has an extremely weak coupling between Cu02 planes, 
and is thus a more exaggerated quasi-2D system than 
most high T c systems; and (iii) this material seems to 
be extremely difficult to dope away from half filling us- 
ing normal preparation techniques (a summary of this 
behaviour may be found in Ref. ]l||). Further, at present 
our experimental understanding of this material is nearly 
complete. Its structure has been determined by x- 
ray and neutron diffraction pLpl to be body-centered- 
tetragonal (K2N1F4 type), and no transition to the or- 
thorhombic phase is observed down to at least 10 K. This 
compound has a three-dimensional AFM structure and 
a Neel temperature Tjv — 251 ± 5 K. In fact, this com- 
pound is considered to be the best available experimental 
realization of the S = h 2D square lattice Heisenberg an- 
tiferromagnet (2DSLHA) pfl, with nearest neighbour 
exchange J = 125 ± 6 meV. Lastly, resistivity measure- 
ments show that it is strongly insulating 0. 

The (essentially) undopable character of this com- 
pound makes it an ideal candidate for angle-resolved 
photoemission (ARPES) experiments. Such studies allow 
for the measurement of the spectral function A'~~'(k, w) 
(simply abbreviated as A(k,u>) from now on) of a single 
hole (left behind by the creation of a photoelectron) prop- 
agating in a Cu02 plane, and such data has recently be- 
come available [ 5| . The results from this experiment thus 
provide an ideal testing ground for theories that purport 
to mimic the low-energy physics of Cu02 planes; e.g., 
the most common model, the so-called t — J model, de- 
scribes spinless vacancies hopping between neighbouring 
sites in a S = \ 2DSLHA. If this model is indeed a good 
representation of a single carrier moving in this plane, it 
should reproduce the full spectral functions measured in 



Ref. B . The point of this paper is to demonstrate how 
well such models of carrier motion in strongly correlated 
systems represent the spectral properties via a compar- 
ison with the full, measured spectral functions of Ref. 



A. Theoretical Comparisons 

Not surprisingly, the ARPES results of Ref. |5) have 
been the subject of intense theoretical studies. For ex- 
ample, it has been noted that the single-hole t— J model 
cannot reproduce the experimental band structure every- 
where in the Brillouin zone, particularly along the (0, 0) 
to (7r, 0), and the (it, 0) to (0, 7r) directions |j. Instead, at 
the very least hoppings beyond near-neighbour must be 
added to the t — J model ||q— [lOf| . The recognition of the 
importance of further hoppings has also been obtained 
in a spin-density wave treatment of the one-band Hub- 
bard model pT| . Work on the more physically plausible 
three-band model has shown that this model is supe- 
rior in reproducing the experimental band structure in 
all regions of the first Brillouin zone []7 12 RJ. However, 
3-band models may be unnecessarily complicated, and 
thus efforts to find a useful and quantitatively accurate 
one-band model are warranted. 

The above comparisons were based on theoretically 
determined and experimentally inferred dispersion rela- 
tions. Of course, experimentally much more information 
is available, namely, the entire spectral function is known 
at many wave vectors. Further, it would be best if the 
comparison of theory and experiment was based on ex- 
act theoretical results. This is not possible unless one 
treats the strong coupling Hamiltonians numerically. Re- 
cently we have completed the first exact diagonalization 
study on a 32-site square lattice for the t — J model 
pi , a formidable task (computationally speaking) . This 
achievement is significant because exact diagonalization 
results obtained with the smaller 4x4 fully square cluster 
are subject to the oddity that this lattice is equivalent to 
the four dimensional hypercube (in the absence of inter- 
actions longer than nearest neighbour). Other exact re- 
sults have been obtained on smaller, non-square clusters, 
or on clusters that lack important wave vectors (namely, 
those along the AFM Brillouin zone boundary); a sum- 
mary of many of these results is contained in a recent 
review jLq]. And lastly and perhaps most importantly, 
while analytical work based on the self-consistent Born 
approximation for the t — J model has been shown to be 
very reliable, at least as far as the dispersion relations 
are concerned E3 , so far no analogous demonstration of 
the accuracy of this approach to include further-than- 
nearest-neighbour hoppings of the t — J model has been 
published. Thus, we consider our numerical, unbiased, 
32-site exact diagonalization approach to a precise man- 
ner in which these models can be compared to experi- 
ment. 



In the present paper we continue our numerical studies 
on the 32-site square lattice, and present a comparison of 
the full spectral functions for the so-called t — t' — t" — J 
model, a generalization of the simpler t — J model that 
includes further hopping processes, as well as three-site, 
spin-dependent hopping. Various groups have claimed 
that this model accurately fits the dispersion relations 
found in the AFM insulator Sr2Cu02Ci2, and such claims 
are something that we can critique with our numerical 
results. Here, besides continuing these comparisons with 
our numerical work, we propose an energy dependence of 
the self-energy for the propagating spin polarons. 



B. Experimental Analysis 

Angle-resolved photoelectron spectroscopy (ARPES) 
has proven to be an invaluable experimental tool for the 
understanding of the electronic structure of the cuprate 
superconductors and related correlated-electron mate- 
rials. The analysis of the ARPES results depends on 
the emission process being primarily direct transitions 
so that the spectra reflect the spectral function A(k, lo). 
This has allowed for the extensive use of ARPES to de- 
termine energy versus k relations for the highest energy 
band in the cuprates, and particularly for determining 
the Fermi surfaces of the metals. Most comparisons be- 
tween calculations and experiment have been performed 
by extracting a quasiparticle energy at each k from both 
the experimental spectra and the calculated A(k, uj) and 
comparing the dispersions thus derived. However, since 
both experiment and theoretical calculations represent 
the spectral function, it should be more instructive to 
compare the two directly rather than just comparing the 
derived dispersions. A direct comparison of spectral func- 
tions has its own difficulties, as are outlined below, but 
hopefully will allow a more complete evaluation of the 
spectral function derived from the Green's function and 
thus avoid the possible errors involved in determining the 
experimental quasiparticle dispersion. 

For two-dimensional materials such as the cuprates, 
the inherent lineshape of the ARPES spectrum may re- 
veal information about the nature of the elementary ex- 
citations. Many previous ARPES studies have examined 
the lineshape of the spectral functions. (For a review of 
this subject, see Ref. jlq ]; for particularly relevant exper- 
imental work, see ]17[[7) In these studies the assumption 
is made that direct transitions dominate the spectral in- 
tensity, and thus one may model the spectral response 
with a single peak, and it has been proposed that the na- 
ture of this peak should be characteristic of a Fermi liq- 
uid, a marginal Fermi liquid, or some other microscopic 
model. Since these fits have all been to a single peak they 
represent only the coherent part of the spectrum. Most 
numerical calculations of A(k, uj) for the cuprates indi- 
cate that there is a great deal of spectral weight in the 
incoherent parts of the excitations. Therefore, attempt- 



ing to fit the cuprate ARPES spectra with a single peak 
of any shape is inherently flawed. In this work we make 
no attempt to deduce a specific lineshape, or broadening 
function, but merely use a few simple fitting procedures 
as discussed below. The essential point is that for the 
purposes of fitting the experimental data, a correct form 
for A(k, lu) is more important than the detailed lineshape 
of a single peak. 

The planar copper oxides are well suited for analy- 
sis since in most of these materials there is a single 
highest energy occupied state (band) which is well sep- 
arated in energy from any other bands. Experimen- 
tally, Sr 2 CuC>2Cl2 has some particular advantages for the 
present study in addition to those mentioned above. It 
is easy to prepare good surfaces for photocmission since 
the material cleaves in a manner similar to mica. In ad- 
dition, the fact that the oxygen content is fixed to the 
stoichiometric value gives us confidence that the actual 
surface under study truly represents a Cu02 plane with a 
single hole per copper site. Thus, while the lineshape and 
broadening appropriate for a particular peak in A(h, uj) 
may be more difficult to determine than for the conduct- 
ing cuprates, it should be more straight forward to de- 
termine the correct Greens function for comparison. 

Our paper is organized as follows. Firstly, in § Q 
we discuss the extraction of the lineshapes of the spin- 
polaron band from the ARPES data. Then, in § IIJ we 
discuss the various modifications to the t — J model nec- 
essary to improve the success of this model in explaining 
the ARPES results. In § IV we present our numerical re- 
sults and the comparison of our work to the experimental 
spectral functions; conclusions are provided in § \V\ 



Since this background function describes the low energy 
valence band, we allow it to be k-dependent. We deter- 
mine this background function at each k by least-square 
fitting, and subtract it off from the experimental data. 

The corrected spectra have the following characteris- 
tics: (i) the peaks broaden and weaken as k moves away 
from the valence band maximum (VBM) at (§,§), (ii) 
the peaks at their sharpest are broad and asymmetri- 
cal, and (iii) although the quasiparticle dispersion is flat 
along the (0, 0) to (it, 0) direction, the spectral weight is 
a maximum at (§, 0). 
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FIG. 1. The background function for the main valence 
band of the insulator (solid line) superimposed on the ARPES 
data, represented by open circles (o), for k = (-|, -|). The re- 
maining spectral weight is associated with the intrinsic spec- 
tral response, coherent as well as incoherent, of a single hole 
in a 2D AFM. 



II. EXTRACTION OF A(k,w) FROM ARPES 
DATA 

The collection of the relevant data has already been 
discussed in Ref. [pi, and in that paper the manner in 
which the quasiparticle dispersion relation was obtained 
from the ARPES data was presented. Here we focus on 
the quasiparticle lineshapes. 

We propose that the experimental data corresponds to 
such a direct-transition-dominated A(k, ui), plus a back- 
ground which influences the low energy part of the data 
and gives rise to an artificially anisotropic peak shape. 
We need to subtract off this background from the experi- 
mental data before they can be compared to the theoret- 
ical A(k,ui). We have little theoretical guidance for the 
background function, but we find that a single Gaussian 
peak plus a constant fits the low energy tail of the ex- 
perimental data very well. In principle, the background 
can also include a function linear in 10 to account for in- 
elastic scattering, and such a linear function should be 
k-independent. However, the best fit at every k gives 
different linear functions. Since this term is quite small 
in all cases, we have chosen to simply set it to be zero. 



As an example of the application of this procedure, 
in Fig. |l| we show the ARPES intensities and the back- 
ground function that we have employed to obtain the 
quasiparticle A(k.,u>) for k = (§,§)• Clearly, the peak 
arising from the main valence band is very well repre- 
sented, even without the inclusion of the inelastic linear 
term — similar satisfactory agreement is found at other 
wave vectors. Thus, we hope that by following this pre- 
scription we are comparing our theoretical work to the 
intrinsic spectral response of a single hole propagating in 
a Cu02 plane in Sr 2 Cu02Cl2. 



III. THEORETICAL MODELS 

To study the theoretical spectrum A(k, ui) of the planar 
copper oxides, we start with the simple t — J model with 
Hamiltonian, 
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where nn means nearest neighbours [Tq] . It is well known 
that this simple, one-band model cannot reproduce all 
the features of the single hole dispersion relation E(k) of 
Sr 2 Cu0 2 Cl 2 as measured by ARPES |Q. While the 
E(k) from the t — J model agrees with experiment along 
the (0, 0) to (tt, tt) direction, two discrepancies remain: 
(i) along the antiferromagnetic Brillouin zone (ABZ) edge 
((tt, 0) to (0, 7r)), the t — J model is less dispersive; and 
(ii) along the (0, 0) to (tt, 0) direction, the model is much 
more dispersive than the experimental result. 

Various modifications to the t — J model have been 
suggested to make the theoretical prediction in better 
agreement with experiment. One suggestion follows from 
the finding that a three-band model works much bet- 
ter than the one-band model [f7|Jl 3|1 , viz. one can mimic 
many features of the three-band model in the one-band 
model by including farther-than-nearest-neighbour hop- 
ping terms. For example, it has become customary to 
include §1 



n 2 = -t' ]T (ct.c^ + H.c.) 

2nn,cr 



(2) 



3nn.er 



where 2nn and 3nn are the second and third-nearest 
neighbours, respectively. Inclusion of Hi significantly 
lowers the anisotropy of the effective mass tensor found in 
the theoretical E(k) around (§, §) |§, but discrepancy 
(ii) above still remains. However, recently it has been 
suggested that so-called three-site, spin-dependent hop- 
ping terms 111] will make the theoretical E(k) from (0, 0) 
to (tt, 0) almost dispersionless PjlCj, thus eliminating this 
problem. These processes are given by 



,/ 
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(3) 



where 5 and 5' are the unit vectors ±x and ±y. 

We have investigated the effects of these terms (H2 
and H3) on -E'(k) of the t — J model using the same 
method as in Ref. pM. In particular, we calculate the 
quasiparticle dispersion relation E(\c) and quasiparticle 
weight Zk of these models for a 32-site square lattice. 
The parameters we employed are t' = 0.3/J, t" = — 0.2i, 
and J = 0.3t, and our results are tabulated in Table 
H (Note that E(k) as defined in Ref. Q is relative to 
the undoped ground state energy.) As pointed out before 
|L4| , the behavior of Z^ of the t — J model along the ABZ 
edge (from (n,0) to (0, n)) does not agree with experi- 
ment. ARPES data shows that the peak intensity is the 
largest at (■§,■§), whereas Z^ for this wave vector for the 
t — J is the smallest along this direction. From Table | it 
is clear that inclusion of farther than nearest neighbour 
hopping terms eliminates this discrepancy. We also note 
that the three-site hopping terms H3 suppress Z^ at all k 



except at the VBM (§,§). Although we cannot produce 
a quantitative comparison of Z^ with ARPES results, we 
will see (in Section fv\ j that this suppression of Z^ makes 
A(k, u>) agree better with ARPES results. 
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FIG. 2. Dispersion relations for single hole propagating in a 
2D square lattice AFM described by the t — t'—t—J model 
(a) without (Hi + Hz) and (b) with (Hi + H 2 + H 3 ) the 
three-site, spin-dependent hopping processes. We have used 
the following parameters: t = — 0.3t, t" = 0.2t, and J = 0.3t. 
The open circles (o) with error bars are the ARPES results 
from Ref. ||. 

To make a comparison of E(k) with experimental re- 
sults, we should express it in units of eV and relative 
to the experimental Fermi energy. This can be done by 
using the energy scale J — 0.125 eV and the location of 
the ARPES peak at k = (§, §). In Fig. | we plot the 
dispersion relation of the Hi + 7Y 2 + H3 with and without 
the three-site hopping terms, and our numerical results 
are compared with ARPES data H . Consistent with pre- 
vious studies g, we find that adding the farther-than- 
nearest-neighbour hopping terms to the t— J model leads 
to an (essentially) isotropic effective mass tensor near the 
VBM. However, there is still some dispersion along the 
(0, 0) to (tt, 0) direction, in contrast to the featureless 
ARPES results. Previous studies ||,|l0| have suggested 
that including the three-site hopping term 7i 3 can re- 
duce the dispersion along this direction, hence making 
the theoretical prediction closer to the experimental re- 
sult. However, from Fig. g(b) we see that H3 only lowers 
the energy of (^, 0) by a small amount (and thus our ex- 



act, unbiased, numerical results are in disagreement with 
earlier analytical work, as well as with work on smaller 
clusters). Despite the lack of perfect agreement along this 
direction, this is the best comparison to the dispersion re- 
lation inferred from the ARPES data for the 7^1+7^2+7^3 
Hamiltonian that we have found. 



IV. CALCULATION OF A(k,w) AND 
COMPARISON WITH ARPES 

The principal quantity of interest in this paper is the 
electron spectral function A(k, w), and is defined by 



A(k,«)=x;i« 
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where Eq and tpQ are the ground-state energy and wave 
function at half filling, respectively, and E^ 1 and tpn -1 
are the energy and wave function of the nth eigenstate of 
the single-hole problem, respectively. To fit the lineshape 
of the ARPES results we need to broaden the delta peaks 
from the theoretical calculation. The continued fraction 
expansion used in Ref. |Q is equivalent to broadening 
the peaks by Lorentzians, 



A(k, W ) = ^|(^- 1 |c k , CT K 
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(5) 



The broadening factor T is a constant, independent of k 
and w. To determine T we fit A(\l,lj) to the ARPES at 
k = (-|, -|) where the peak is most well defined. For the 
best agreement, we have to use T = 0.6t ~ 250 meV. 
Note that this is as large as the coherent bandwidth, 
which is 280 ± 60 meV. Now, since all energy scales and 
r are fixed by the spectrum at k = (^, |), no fitting is 
needed for other k. 

In Fig. g we plot A(k, uj) on top of the corrected 
ARPES data at the available k from (0,0) to (tt,tt). 
(These experimental data are probably the most reliable 
because in the experiment the angle between the elec- 
tric field vector, the sample axes, and the ejected elec- 
trons remains constant.) Since the theoretical model is 
supposed to describe only the low energy physics of the 
CuC>2 plane, we truncate the spectrum at W = — 1.8 eV. 
It is obvious that the model can reproduce the experi- 
mental finding that the spectral weight is maximum at 
the VBM and that the peak broadens and weakens as k 
moves away from this point in the Brillouin zone. Fur- 
ther, we find that the incoherent part of the theoretical 
spectra make the peaks asymmetric, also in qualitative 
agreement with experiment. However, they are not as 
asymmetric as the experimental results, and the compar- 
ison is not entirely encouraging. Of course, the question 
is: is this the best this theoretical model can do? 
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FIG. 3. Comparison of the theoretical A(k,ui) for 
TCi + Ti.2 + Ti-3 to the ARPES data at available wavevectors k 
from (0,0) to (tt, 7r). Solid lines are A(k,u) calculated using 
constant broadening F = 0.6t (refer to Eq. (H)) at the k in- 
dicated in the graph. Open circles (o) are corrected ARPES 
data at the same k. (Since the experiments were not per- 
formed at exactly the same wave vectors that we can treat nu- 
merically, for the theoretical k = (f, f) we show the ARPES 
data at (0.37T, 0.37r), and for the theoretical k = (2-p, 2p) we 
show the ARPES data at (0.8tt, 0.8tt).) 

Looking at the VBM peak, which is the sharpest and 
the most prominent, it is clear that we cannot completely 
fit the ARPES result with simple Lorentzians, since the 
high energy part of the peak falls off much faster than a 
Lorentzian. Contrary to the ARPES of metals, this faster 
fall off cannot be accounted for by including a Fermi func- 
tion. To improve the comparison of our numerics and 
experiment, without microscopic justification we use an 
energy-dependent broadening function T(u>). Analogous 
to Eq. (H), we write 



A(k, W )=^|(^- 1 |c k , CT K)| i 
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where Af(E) is the normalization factor, 

I» 
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From the spectrum at (?, |) in Fig. |[ the low energy tail 
(incoherent part) of the spectrum seems to need larger 
broadening. For simplicity we choose T(uj) to be linear 
in lo: 



r(w) = max(0, -0.176 eV - 0.55w) 



(8) 



Note that at the first peak {w 0.85 eV), T(u>) ~ 290 

meV which is similar to the constant T mentioned above. 



along the (0,0) to (tt, 0), and the (tt, 0) to (0, n) direc- 
tions, respectively. Figures || ||, and || show that our 
model fits the ARPES lineshape at many k very well. 
The least satisfactory fit is at k = (§,0), which is also 
the least satisfactory k point in the dispersion relation 
comparison, as seen in Fig. || Despite this, our model 
does produce the largest spectral weight (c./. Table |) 
along the (0, 0) to (ir, 0) direction at this wave vector, 
similar to experiment. 
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FIG. 4. This figure shows the comparison of our spec- 
tral functions to wave vectors along (k,k), the same 
as in Fig. pi except that A(k,u>) is calculated using 
r(w) = max(0, —0.176 eV — 0.55a;). We have also included 
solid arrows to indicate the positions of the first peaks of the 
ARPES data as shown in Fig. H, and hollow arrows to indicate 
the numerically determined E(k). 

In Fig. [I] we plot A(k, uj) calculated from Eqs. (g) and 
(0) on top of the corrected ARPES data for those k from 
(0,0) to (7r,7r). This shows that our energy-dependent 
broadening function is able to produce a much improved 
fit at all k. In Figs. [5] and || we plot the comparison 
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FIG. 5. Same as Fig. | but along the (0, 0) to (tt, 0) direc- 
tion. 

In principle, we can allow the broadening function to 
be dependent on k, r(k, u>) [0. For example, instead 
of determining the two coefficients in Eq.(0) using the 
spectrum at (f,5), we can determine them at each k. 
We have completed such fittings, but we do not find much 
better results, and in particular, the discrepancy between 
theoretical and experimental spectra at (§, 0) remains. 



V. CONCLUSIONS 

To summarize, to compare the theoretical and experi- 
mental lineshapes we have argued that the ARPES data 
of Sr 2 Cu02Cl2 is dominated by direct transitions, and 



thus with the valence band contribution appropriately 
subtracted off, one can obtain the spectral function for a 
single hole propagating in a CuC>2 plane. From these ex- 
perimental spectral functions it is obvious that the spec- 
trum at each k cannot be fit by a single, narrow peak; 
instead, the incoherent part is important. 
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not be seen as direct support for the marginal Fermi liq- 
uid hypothesis, since we are studying an insulator and 
thus we cannot predict that this linear term survives in 
the (anomalous) metallic state. 

The least satisfactory fit is found along the (ir, 0) to 
(0,0) direction, and at (§,0) our calculated spectrum is 
significantly lower than the ARPES result. To date, only 
three-band models have provided a reasonable agreement 
with the dispersion near this point. 

An interesting question that follows from our results 
is: What is the origin of the incoherent processes that 
contribute so strongly to spectral functions of a sin- 
gle hole propagating in a strongly correlated, half-filled 
state? Two possibilities seem worthy of further study. 
One, long-wavelength spin excitations, excluded from 
our work due to our use of a finite cluster, are impor- 
tant. Two, some other excitations, e.g. lattice waves, 
are strongly coupled to either the hole motion or the 
spin waves, or possibly both, and provide an effective 
energy-dependent damping of these quasiparticles. 
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FIG. 6. Same as Fig. |j but along the (n, 0) to (0, tv) direc- 
tion. (Similar to the circumstances discussed in Fig. fa, for 



the theoretical k = (^ L , f ) we show the ARPES results for 
k = (0.37T, 0.771"), the latter found from an average of the data 
for (0.7tt, 0.3tt) and (0.3tt, 0.7tt).) 



We have studied the effects of including farther-than- 
nearest neighbour and three-site, spin-dependent hop- 
ping terms to the t — J model. Consistent with previous 
studies, we find that these terms increase the dispersion 
along the antiferromagnetic Brillouin zone edge ((77, 0) to 
(0, 7r)), and suppress that along the (n, 0) to (0, 0) direc- 
tion. In addition, these terms suppress the quasi-particle 
weight except at the valence band maximum (§j £)■ This 
is shown to be essential in fitting the theoretical lineshape 
to the experimental results. 

We also find that if we use an energy-dependent broad- 
ening function T(lo), we can obtain a reasonably good fit 
in an energy range of up to 1.8 eV below the Fermi level. 
The fact that we find a damping linear in energy should 
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TABLE I. Quasiparticle weight Z^ and energy E(k) of the 
t-t'-t"-J model with {Hi+Hi+H-j,) and without (H1+H2) 
the three-site, spin-dependent hopping terms, calculated on 
a 32-site square lattice. The parameters that we used are 
t' = -0.3£, t" = 0.2t, and J = 0.3t E(k) is in units of t 
and is measured relative to the half-filled ground state energy 
Eo = -11.3297£. 
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